Manipulating the diffusion energy barrier at the lithium metal electrolyte interface for dendrite-free long-life batteries.

Constructing an artificial solid electrolyte interphase (SEI) on lithium metal electrodes is a promising approach to address the rampant growth of dangerous lithium morphologies (dendritic and dead Li0) and low Coulombic efficiency that plague development of lithium metal batteries, but how Li+ transport behavior in the SEI is coupled with mechanical properties remains unknown. We demonstrate here a facile and scalable solution-processed approach to form a Li3N-rich SEI with a phase-pure crystalline structure that minimizes the diffusion energy barrier of Li+ across the SEI. Compared with a polycrystalline Li3N SEI obtained from conventional practice, the phase-pure/single crystalline Li3N-rich SEI constitutes an interphase of high mechanical strength and low Li+ diffusion barrier. We elucidate the correlation among Li+ transference number, diffusion behavior, concentration gradient, and the stability of the lithium metal electrode by integrating phase field simulations with experiments. We demonstrate improved reversibility and charge/discharge cycling behaviors for both symmetric cells and full lithium-metal batteries constructed with this Li3N-rich SEI. These studies may cast new insight into the design and engineering of an ideal artificial SEI for stable and high-performance lithium metal batteries.

Detonation Performance of Insensitive Nitrogen-Rich Nitroenamine Energetic Materials Predicted from First-Principles Reactive Molecular Dynamics Simulations.

Because of the excellent combination of high detonation and low sensitivity properties of the 1,1-diamino-2,2-dinitroethylene (FOX-7) energetic material (EM), it is useful to explore new energetic derivatives that start with the FOX-7 structure. However, most such derivatives are highly sensitive, making them unsuitable for EM applications. An exception is the new nitroenamine EM, 1,1-diamino-2-tetrazole-2-nitroethene (FOX-7-T) (synthesized by replacing a nitro group with a tetrazole ring), which exhibits good stability. Unfortunately, FOX-7-T shows an unexpected much lower detonation performance than FOX-7, despite its higher nitrogen content. To achieve an atomistic understanding of the insensitivity and detonation performance of FOX-7 and FOX-7-T, we carried out reactive molecular dynamics (RxMD) using the ReaxFF reactive force field and combined quantum mechanics MD (QM-MD). We found that the functional group plays a significant role in the initial decomposition reaction. For FOX-7, the initial decomposition involves only simple hydrogen transfer reactions at high temperature, whereas for FOX-7-T, the initial reaction begins at much lower temperature with a tetrazole ring breaking to form N2, followed by many subsequent reactions. Our first-principles-based simulations predicted that FOX-7-T has 34% lower CJ pressure, 15% lower detonation velocity, and 45% lower CJ temperature than FOX-7. This is partly because a larger portion of the FOX-7-T mass gets trapped into condensed phase carbon clusters at the CJ point, suppressing generation of gaseous CO2 and N2 final products, leading to reduced energy delivery. Our findings suggest that the oxygen balance is an important factor to be considered in the design of the next generation of high-nitrogen-containing EMs.

Molecular Dynamics Studies of Atomistically Determined Fibrillar Assemblies: Comparison of the Rippled ß-Sheet, Pleated ß-Sheet, and Herringbone Structures.

Pauling and Corey expected that a racemic mixture would result in a rippled ß-sheet, however, it has been known from experiments that the racemic mixtures of triphenylalanine lead to a herringbone structure. Because of the theoretical limitations concerning crystal structures such as rippled ß-sheet, it is inevitable to understand how the interplay of the amino acids prefers a specific structural motif. In this paper we use molecular dynamics to understand the sequence- and enantiomer-dependent structures by comparisons between rippled ß-sheet and pleated ß-sheet, solvated and anhydrous rippled ß-sheet, and rippled ß-sheet and the herringbone structure, based on thermodynamics and structures at the atomic level. The tripeptides select the favored structure that can be stabilized through aromatic or hydrogen bonding interactions between tripeptides. Furthermore, the solubility is determined by the environment of space that is created around the side chains. Our findings provide comprehensive insight into the crystallized fibril motif of the polypeptide.

Argentophilic Interactions, Flexibility, and Dynamics of Pyrrole Cages Encapsulating Silver(I) Clusters.

Recently, pyrrole cages have been synthesized that encapsulate ion pairs and silver(I) clusters to form intricate supramolecular capsules. We report here a computational analysis of these structures using density functional theory combined with a semiempirical tight-binding approach. We find that for neutral pyrrole cages, the Gibbs free energies of formation provide reliable predictions for the ratio of bound ions. For charged pyrrole cages, we find strong argentophilic interactions between Ag ions on the basis of the calculated bond indices and molecular orbitals. For the cage with the Ag4 cluster, we find two minimum-geometry conformations that differ by only 6.5 kcal/mol, with an energy barrier <1 kcal/mol, suggesting a very flexible structure as indicated by molecular dynamics. The predicted energies of formation of [Agn⊂1]n-3+ (n = 1-5) cryptands provide low energy barriers of formation of 5-20 kcal/mol for all cases, which is consistent with the experimental data. Furthermore, we also examined the structural variability of mixed-valence silver clusters to test whether additional geometrical conformations inside the organic cage are thermodynamically accessible. In this context, we show that the time-dependent density functional theory UV-vis spectra may potentially serve as a diagnostic probe to characterize mixed-valence and geometrical configurations of silver clusters encapsulated into cryptands.

Enhancing Multifunctionality: Optimal Properties of Iron-Oxide-Reinforced Polyvinylidene Difluoride Unveiled Through Full Atom Molecular Dynamics Simulations.

Nanocomposites made of magnetite (Fe3O4) nanoparticles (NP)s with different surface chemistry and polyvinyl difluoride (PVDF) polymer were investigated using full atom molecular dynamics (MD) simulation. NPs with hydroxyl (OH), hexanoic, and oleic acid terminations were considered in this study. The effect of each surface chemistry was investigated in terms of the mechanical properties, the distribution of the internal energy around the NP, and the chain polarization gradient from the interface to the bulk. From this investigation, we find that oleic acid termination, although the most popular, is less favorable for interfacial interaction and local polarization. The OH-terminated NP results in the best configuration for the properties investigated. The hexanoic acid-grafted NP presents a good compromise. Hydrogen bonding governs the induced response of the nanocomposites. Although the hexanoic acid grafted NP presents less hydrogen bonding than the OH-terminated case, the conformation of the hexanoic acid acts as a mobility flow inhibitor, leading to a performance comparable to that of the OH-terminated NP composite. This work led to investigating routes to make nanocomposite materials with optimized properties. These results shed light on the multiple combinations offered by nanocomposites that go beyond the conventional effects of size.

Facet-Dependent Oxygen Evolution Reaction Activity of IrO2 from Quantum Mechanics and Experiments.

The diversity of chemical environments present on unique crystallographic facets can drive dramatic differences in catalytic activity and the reaction mechanism. By coupling experimental investigations of five different IrO2 facets and theory, we characterize the detailed elemental steps of the surface redox processes and the rate-limiting processes for the oxygen evolution reaction (OER). The predicted complex evolution of surface adsorbates and the associated charge transfer as a function of applied potential matches well with the distinct redox features observed experimentally for the five facets. Our microkinetic model from grand canonical quantum mechanics (GC-QM) calculations demonstrates mechanistic differences between nucleophilic attack and O-O coupling across facets, providing the rates as a function of applied potential. These GC-QM calculations explain the higher OER activity observed on the (100), (001), and (110) facets and the lower activity observed for the (101) and (111) facets. This combined study with theory and experiment brings new insights into the structural features that either promote or hinder the OER activity of IrO2, which are expected to provide parallels in structural effects on other oxide surfaces.

CO2-Promoted Electrocatalytic Reduction of Chlorinated Hydrocarbons.

Electrochemical reactions and their catalysis are important for energy and environmental applications, such as carbon neutralization and water purification. However, the synergy in electrocatalysis between CO2 utilization and wastewater treatment has not been explored. In this study, we find that the electrochemical reduction of chlorinated organic compounds such as 1,2-dichloroethane, trichloroethylene, and tetrachloroethylene into ethylene in aqueous media, which is a category of challenging reactions due to the competition of H2 evolution, can be substantially enhanced by simultaneously carrying out the reduction of CO2 on an easily prepared and cost-effective Cu metal catalyst. In the case of 1,2-dichloroethane dechlorination, a 6-fold improvement in Faradaic efficiency and a 19-fold increase in partial current density are demonstrated. Through electrochemical kinetic studies, in situ Raman spectroscopy, and computational simulations, we further find that CO2 reduction reduces hydrogen coverage on the Cu catalyst, which not only exposes more active sites for the dechlorination reaction but also enhances the effective reductive potential on the catalyst surface and reduces the kinetic barrier of the rate-determining step.

Methotrexate Inhibits the Binding of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Receptor Binding Domain to the Host-Cell Angiotensin-Converting Enzyme-2 (ACE-2) Receptor.

As the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus mutates, finding effective drugs becomes more challenging. In this study, we use ultrasensitive frequency locked microtoroid optical resonators in combination with in silico screening to search for COVID-19 drugs that can stop the virus from attaching to the human angiotensin-converting enzyme 2 (hACE2) receptor in the lungs. We found 29 promising candidates that could block the binding site and selected four of them that were likely to bind very strongly. We tested three of these candidates using frequency locked optical whispering evanescent resonator (FLOWER), a label-free sensing method based on microtoroid resonators. FLOWER has previously been used for sensing single macromolecules. Here we show, for the first time, that FLOWER can provide accurate binding affinities and sense the inhibition effect of small molecule drug candidates without labels, which can be prohibitive in drug discovery. One of the candidates, methotrexate, showed binding to the spike protein 1.8 million times greater than that to the receptor binding domain (RBD) binding to hACE2, making it difficult for the virus to enter cells. We tested methotrexate against different variants of the SARS-CoV-2 virus and found that it is effective against all four of the tested variants. People taking methotrexate for other conditions have also shown protection against the original SARS-CoV-2 virus. Normally, it is assumed that methotrexate inhibits the replication and release of the virus. However, our findings suggest that it may also block the virus from entering cells. These studies additionally demonstrate the possibility of extracting candidate ligands from large databases, followed by direct receptor-ligand binding experiments on the best candidates using microtoroid resonators, thus creating a workflow that enables the rapid discovery of new drug candidates for a variety of applications.

Electrochemical Nitrate Reduction Catalyzed by Two-Dimensional Transition Metal Borides.

We investigated 2D transition metal borides (MBenes) for the efficient conversion of nitrate to ammonia. MBenes have been previously shown to bind oxygen with extraordinary strength, which should translate toward selective adsorption of nitrate in aqueous media. Using Density Functional Theory, we screened MBenes by computing their nitrate and water adsorption energies, seeking materials with strong nitrate binding and weak water binding. We identified MnB, CrB, and VB as the best materials for selective nitrate adsorption and proceeded by computing their free energies for generating ammonia. Of the three candidates, CrB requires the lowest overpotential, making it the best candidate. To further decrease the overpotential, we doped the CrB MBene with secondary transition metals and found the addition of Mn to the active site further reduced the overpotential. We then computed the reaction mechanism grand canonically to observe the effect of applied potential on the free energy landscape.

How Permanent Are the Permanent Macrodipoles of Anthranilamide Bioinspired Molecular Electrets?

Dipoles are ubiquitous, and their impacts on materials and interfaces affect many aspects of daily life. Despite their importance, dipoles remain underutilized, often because of insufficient knowledge about the structures producing them. As electrostatic analogues of magnets, electrets possess ordered electric dipoles. Here, we characterize the structural dynamics of bioinspired electret oligomers based on anthranilamide motifs. We report dynamics simulations, employing a force field that allows dynamic polarization, in a variety of solvents. The results show a linear increase in macrodipoles with oligomer length that strongly depends on solvent polarity and hydrogen-bonding (HB) propensity, as well as on the anthranilamide side chains. An increase in solvent polarity increases the dipole moments of the electret structures while decreasing the dipole effects on the moieties outside the solvation cavities. The former is due to enhancement of the Onsager reaction field and the latter to screening of the dipole-generated fields. Solvent dynamics hugely contributes to the fluctuations and magnitude of the electret dipoles. HB with the solvent weakens electret macrodipoles without breaking the intramolecular HB that maintains their extended conformation. This study provides design principles for developing a new class of organic materials with controllable electronic properties. An animated version of the TOC graphic showing a sequence of the MD trajectories of short and long molecular electrets in three solvents with different polarities is available in the HTML version of this paper.

2D Metal/Graphene and 2D Metal/Graphene/Metal Systems for Electrocatalytic Conversion of CO2 to Formic Acid.

Efficiently transforming CO2 into renewable energy sources is crucial for decarbonization efforts. Formic acid (HCOOH) holds great promise as a hydrogen storage compound due to its high hydrogen density, non-toxicity, and stability under ambient conditions. However, the electrochemical reduction of CO2 (CO2 RR) on conventional carbon black-supported metal catalysts faces challenges such as low stability through dissolution and agglomeration, as well as suffering from high overpotentials and the necessity to overcome the competitive hydrogen evolution reaction (HER). In this study, we modify the physical/chemical properties of metal surfaces by depositing metal monolayers on graphene (M/G) to create highly active and stable electrocatalysts. Strong covalent bonding between graphene and metal is induced by the hybridization of sp and d orbitals, especially the sharp d z 2 ${{d}_{{z}^{2}}}$ , d y z ${{d}_{yz}}$ , and d x z ${{d}_{xz}}$ orbitals of metals near the Fermi level, playing a decisive role. Moreover, charge polarization on graphene in M/G enables the deposition of another thin metallic film, forming metal/graphene/metal (M/G/M) structures. Finally, evaluating overpotentials required for CO2 reduction to HCOOH, CO, and HER, we find that Pd/G, Pt/G/Ag, and Pt/G/Au exhibit excellent activity and selectivity toward HCOOH production. Our novel 2D hybrid catalyst design methodology may offer insights into enhanced electrochemical reactions through the electronic mixing of metal and other p-block elements.

Atomistic mechanisms underlying plastic flow at ultralow yield stress in ductile carbon aerogels.

We investigated carbon aerogel samples with super low densities of 0.013 g cm-3 (graphite is 2.5) and conducted compression experiments showing a very low yield stress of 5-8 kPa. To understand the atomistic mechanisms operating in these super low density aerogels, we present a computational study of the mechanical response of very low-density amorphous carbonaceous materials. We start from our previously derived atomistic models (based on the DynReaxMas method) with a density of 0.16 g cm-3 representing the core regions of carbon aerogels. We considered three different phases exhibiting either a fiber-like clump morphology interconnected with string-like units or a more reticulated framework. We subjected these phases to compression and shear deformations and analyzed the resulting plastic response via an inherent-structure protocol. Strikingly, we find that these materials possess shear plastic relaxation modes with extremely low values of yield stress, negligible with respect to the finite values predicted outside this "zero-stress" region. This is followed by a succession of two additional regimes with increasing yield stress values. Our analysis of the atomistic relaxation mechanisms finds that these modes have a collective and cooperative character, taking the form of nanoscopic shear bands within the clumps. These findings rationalize our experimental observations of very low-stress plastic deformation modes in carbon aerogels, providing the first steps for developing a predictive multi-scale modeling of the mechanical properties of aerogel materials.

Understanding Ionic Diffusion Mechanisms in Li2S Coatings for Solid-State Batteries: Development of a Tailored Reactive Force Field for Multiscale Simulations.

In order to investigate Li2S as a potential protective coating for lithium anode batteries using superionic electrolytes, we need to describe reactions and transport for systems at scales of >10,000 atoms for time scales beyond nanoseconds, which is most impractical for quantum mechanics (QM) calculations. To overcome this issue, here, we first report the development of the reactive analytical force field (ReaxFF) based on density functional theory (DFT) calculations on model systems at the PBE0/TZVP and M062X/TZVP levels. Then, we carry out reactive molecular dynamics simulations (RMD) for up to 20 ns to investigate the diffusion mechanisms in bulk Li2S as a function of vacancy density, determining the activation barrier for diffusion and conductivity. We show that RMD predictions for diffusion and conductivity are comparable to experiments, while results on model systems are consistent with and validated by short (10-100 ps) ab initio molecular dynamics (AIMD). This new ReaxFF for Li2S systems enables practical RMD on spatial scales of 10-100 nm (10,000 to 10 million atoms) for the time scales of 20 ns required to investigate predictively the interfaces between electrodes and electrolytes, electrodes and coatings, and coatings and electrolytes during the charging and discharging processes.

Origins of Enhanced Enantioselectivity in the Pd-Catalyzed Decarboxylative Allylic Alkylation of N-Benzoyl Lactams.

We explore the origins of the marked improvement in enantioselectivity in the inner-sphere (PHOX)Pd-catalyzed allylic alkylation of N-benzoyl lactam nucleophiles over their carbocyclic counterparts. We employ density functional theory calculations to aid in the interpretation of experimental results. Ultimately, we propose that the enhancement in enantioselectivity arises primarily from noncovalent interactions between the substrate and ligand rather than secondary substrate chelation, as previously hypothesized.

Fully activated structure of the sterol-bound Smoothened GPCR-Gi protein complex.

Smoothened (SMO) is an oncoprotein and signal transducer in the Hedgehog signaling pathway that regulates cellular differentiation and embryogenesis. As a member of the Frizzled (Class F) family of G protein-coupled receptors (GPCRs), SMO biochemically and functionally interacts with Gi family proteins. However, key molecular features of fully activated, G protein-coupled SMO remain elusive. We present the atomistic structure of activated human SMO complexed with the heterotrimeric Gi protein and two sterol ligands, equilibrated at 310 K in a full lipid bilayer at physiological salt concentration and pH. In contrast to previous experimental structures, our equilibrated SMO complex exhibits complete breaking of the pi-cation interaction between R4516.32 and W5357.55, a hallmark of Class F receptor activation. The Gi protein couples to SMO at seven strong anchor points similar to those in Class A GPCRs: intracellular loop 1, intracellular loop 2, transmembrane helix 6, and helix 8. On the path to full activation, we find that the extracellular cysteine-rich domain (CRD) undergoes a dramatic tilt, following a trajectory suggested by positions of the CRD in active and inactive experimental SMO structures. Strikingly, a sterol ligand bound to a shallow transmembrane domain (TMD) site in the initial structure migrates to a deep TMD pocket found exclusively in activator-bound SMO complexes. Thus, our results indicate that SMO interacts with Gi prior to full activation to break the molecular lock, form anchors with Gi subunits, tilt the CRD, and facilitate migration of a sterol ligand in the TMD to an activated position.

Heck Migratory Insertion Catalyzed by a Single Pt Atom Site.

Single-atom catalysts (SACs) have generated excitement for their potential to downsize metal particles to the atomic limit with engineerable local environments and improved catalytic reactivities and selectivities. However, successes have been limited to small-molecule transformations with little progress toward targeting complex-building reactions, such as metal-catalyzed cross-coupling. Using a supercritical carbon-dioxide-assisted protocol, we report a heterogeneous single-atom Pt-catalyzed Heck reaction, which provides the first C-C bond-forming migratory insertion on SACs. Our quantum mechanical computations establish the reaction mechanism to involve a novel C-rich coordination site (i.e., PtC4) that demonstrates an unexpected base effect. Notably, the base was found to transiently modulate the coordination environment to allow migratory insertion into an M-C species, a process with a high steric impediment with no previous example on SACs. The studies showcase how SACs can introduce coordination structures that have remained underexplored in catalyst design. These findings offer immense potential for transferring the vast and highly versatile reaction manifold of migratory-insertion-based bond-forming protocols to heterogeneous SACs.

Quantum mechanics based non-bonded force field functions for use in molecular dynamics simulations of materials and systems: The nitrogen and oxygen columns.

Accurate Force Fields (FFs) are essential for Molecular Dynamics (MD) simulations of the dynamics of realistic materials in terms of atomic-level interactions. The FF parameters of short-range valence interactions can be derived through Quantum Mechanical (QM) calculations on model systems practical for QM (<300 atoms). Similarly, the dynamic electrostatic interactions can be described with methods such as QEq or PQEq that allow charges and polarization to adjust dynamically. However, accurately extracting long-range van der Waals (vdW) interactions from QM calculations poses challenges due to the absence of a definitive method to distinguish between the different energetic components of electrostatics, polarization, vdW, hydrogen bonding, and valence interactions. To do this we use the Perdew-Burke-Ernzerhof flavor of Density Functional Theory, including empirical D3 vdW corrections, to predict the Equation of State for each element (keeping any covalent bonds fixed), from which we obtain the two-body vdW nonbond potential. Here, we extend these calculations to include non-bonded parameters for the N and O columns of the periodic table so that we now describe columns 15 (N), 16 (O), 17 (F), and 18 (Ne) of the periodic table. For these 20 elements, we find that the two-body vdW potentials can all be mapped to a single universal two-body curve, with just three scaling parameters: Re, De, and L. We refer to this as the Universal NonBond (UNB) potential. We expect this to be useful for new MD simulations and a helpful starting point to obtain UNB parameters for the remainder of the periodic table.

Genipin and pyrogallol: Two natural small molecules targeting the modulation of disordered proteins in Alzheimer's disease.

Alzheimer's disease (AD) is characterized by the aggregation of disordered proteins, such as amyloid beta (Aß) and tau, leading to neurotoxicity and disease progression. Despite numerous efforts, effective inhibitors of Aß and tau aggregates have not been developed. Thus, we aimed to screen natural small molecules from crude extracts that target various pathologies and are prescribed for patients with neurological diseases. In this study, we screened 162 natural small molecules prescribed for neurological diseases and identified genipin and pyrogallol as hit compounds capable of simultaneously regulating the aggregation of Aß and tau K18. Moreover, we confirmed the dual modulatory effects of these compounds on the reduction of amyloid-mediated neurotoxicity in vitro and the disassembly of preformed Aß42 and tau K18 fibrils. Furthermore, we observed the alleviatory effects of genipin and pyrogallol against AD-related pathologies in triple transgenic AD mice. Molecular dynamics and docking simulations revealed the molecular interaction dynamics of genipin and pyrogallol with Aß42 and tau K18, providing insights into their suppression of aggregation. Our findings suggest the therapeutic potential of genipin and pyrogallol as dual modulators for the treatment of AD by inhibiting aggregation or promoting dissociation of Aß and tau.

Discovery of 3-Phenyl Indazole-Based Novel Chemokine-like Receptor 1 Antagonists for the Treatment of Psoriasis.

Chemokine-like receptor 1 (CMKLR1)─a G protein-coupled receptor─has functional roles in the immune system and related diseases, including psoriasis and metabolic diseases. Psoriasis is a chronic inflammatory disease characterized by skin redness, scaliness, and itching. In this study, we sought to develop novel CMKLR1 antagonists by screening our in-house GPCR-targeting compound library. Moreover, we optimized a phenylindazole-based hit compound with antagonistic activities and evaluated its oral pharmacokinetic properties in a murine model. A structure-based design on the human CMKLR1 homology model identified S-26d as an optimized compound that serves as a potent and orally available antagonist with a pIC50 value of 7.44 in hCMKLR1-transfected CHO cells. Furthermore, in the imiquimod-induced psoriasis-like mouse model, oral administration of S-26d for 1 week significantly alleviated modified psoriasis area and severity index scores (severity of erythema, scaliness, skin thickness) compared with the control group.

Functional Group-Dependent Proton Conductivity of Phosphoric Acid-Doped Ion-Pair Coordinated Polymer Electrolytes: A Molecular Dynamics Study.

Toward deployment of high-temperature polymer electrolyte membrane fuel cells (HT-PEMFCs) in our daily lives, multiple research efforts have been dedicated to develop high-performance phosphate-doped polymer electrolytes. Recently, ion-pair coordinated polymers have garnered attention for their high stability and proton conductivity. However, a comprehensive understanding of how proton transport properties are modified by the functional groups present in these polymers is still lacking. In this study, we employ molecular dynamics (MD) simulations to investigate the impact of different functional group types and conversion ratios on conductivity. We find that Grotthuss-type hopping transport predominantly governs the overall conductivity, surpassing vehicular transport by factors of 100-1000. As conductivity scales with proton concentration, we observe that less-bulky functional groups offer advantages by minimizing the volume expansion associated with increased conversion ratios. Additionally, we show that a strong ion-pair interaction between the cationic functional group and the phosphate anion disrupts the suitable intermolecular orientations required for efficient proton hopping between phosphate and phosphoric acid molecules, thereby diminishing the proton conductivity. Our study underscores the importance of optimizing the strength of ion-pair interactions to balance stability and proton conductivity, thus paving the way for the development of ion-pair coordinated polymer electrolytes with improved performance.